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Abstract 

Active cancellation of broadband random noise requires the detection of the incoming noise 
with some time advance. In an duct for example this advance must be larger than the delays in 
the secondary path from the control source to the error sensor In this paper it is shown that, in 
some cases, the advance required for perfect noise cancellation is theoretically infinite because 
the inverse of the secondary path, which is required for control, can include an infinite non- 
causal response. This is shown to be the result of two mechanisms: in the single-channel case 
(one control source and one error sensor), this can arise because of strong echoes in the control 
path. In the multi-channel case this can arise even in free field simply because of an unfortunate 
placing of sensors and actuators. In the present paper optimal feedforward control is derived 
through analytical and numerical computations, in the time and frequency domains. It is shown 
that, in practice, the advance required for significant noise attenuation can be much larger than the 
secondary path delays. Practical rules are also suggested in order to prevent infinite non-causaUty 
from appearing. 



1 Introduction 

Active Noise Control (ANC) can help in reducing broadband random noises such as road noise inside 
cars (see e.g. Refs. [1-3]), fan noise travelling in ducts (Ref. IH) or noise induced by the boundary- 
layer pressure fluctuations around aircraft fuselages (Refs. [5-6]). However, in the broadband random 
case, the noise arriving at the minimization microphones is mostly unpredictable, and the constraint 
of causality, which applies to the control filters, can resti^ain the performances of ANC (cf. Ref. |7 1). 
Therefore, in practice, ANC set-ups often rely on feedforward control which involves one or several 
reference sensors detecting the incoming noise before it reaches the minimization microphones. In 
the single-channel case of noise travelling along one direction, the constraint of causality is easy to 
interpret in terms of acoustic propagation: the primary noise must be detected before it reaches the 
minimization sensor, with a time advance larger than the propagation time from the control source to 
the minimization microphone. The time advance must also compensate the electrical and computation 
delays in the control processing unit. However, in more complex cases, it is not easy to link the 
causality constraint to propagation delays. Indeed, no design rule guarantees the causality of perfect 
noise cancellation without effective computation of the optimal control filters. 
In this context, this paper shows that perfect noise cancellation can require control filters with an 
infinite non-causal response, which means that perfect feedforward control can require the detection 
of the incoming noise with an infinite time advance. Analytical computations on simple ANC set-ups 
show that this infinite non-causality is the theoretical result of two mechanisms: in the single-channel 
case of one control source and one minimization microphone, it is the consequence of strong echoes in 
the secondary path. In the multi-channel case, it can arise even in free field because of an unfortunate 
geometrical arrangement of the actuators and sensors. In the paper, the performance of causally- 
constrained ANC is also investigated through numerical computations. It is shown that, in practice. 



achieving a significant noise attenuation can require the detection of the incoming noise with a very 
large time advance. 

The study presented in this paper was motivated by the unexpected difficulties encountered when 
implementing multi-channel broadband ANC in a large room at LMA. Although many successful 
experiments of random noise control in enclosures had been reported (see e.g. Refs. fSlj and f9\), it 
appeared that efficient control in the LMA room required feeding the reference signal to the control 
unit with a time advance of several seconds (see Ref. [10|). Previous papers have reported that the 
response of optimal control filters could include a very long non-causal component (e.g. in Refs. 11 11. 
fT2l and |[T3l ). but, although a link between non-minimum phase responses and reverberation had 
been observed earlier (see Ref. (T4i ). this has been seen in the context of ANC as a purely algebraic 
problem. The present paper is an attempt to give an acoustic interpretation of control non-causality 
in order to derive design guidelines for ANC set-ups. 

In section |2lof this paper, the causality of ANC is studied in the single-channel case. After recalling 
the case of noise cancellation in an infinite duct, it is shown that one echo in the secondary path can 
lead to infinitely non-causal control. Optimal non-causal control and causally-constrained control 
are then derived for the case of a single monopole source and a single pressure microphone in the 
corner of two rigid walls. Numerical simulations in the frequency and time domain show that, in 
some cases, the incoming primary noise must be detected with a very large time advance for efficient 
control. Finally it is shown that a slight change in the actuator and sensor locations can relax the 
causality constraint, which suggests guidelines for the design of ANC set-ups not prone to infinite 
non-causality. 

In section |3j ANC causality is investigated in the multi-channel case. Optimal non-causal con- 
trol is derived analytically for a free-field set-up including two monopole sources and two micro- 
phones, and some arrangements are shown to lead to infinitely non-causal noise cancellation. The 



performance of causally-constrained control is evaluated through numerical simulations involving 
a Filtered-Reference Least Mean Square (FXLMS) algorithm. In this 2x2 case it is shown that 
non-causality is critical when the condition number of the matrix of secondary paths is high at low- 
frequency. 

The inversion of acoustic paths, which is required for real-time ANC, is also needed in devices for 
"virtual acoustics" which aim at a sound reproduction conveying accurate information such as the 
location of the recorded noise sources (see Ref. 115.1 ). Since sound reproduction is not a real-time 
process, the non-causality arising from the acoustic path inversion is usually not critical. However 
this paper shows that care must be taken for this inversion, even in the context of virtual acoustics, 
because the additional advance required for causal inversion can be much larger than the delays in the 
direct path. 

2 Non-causality in the single-channel case with reflections 
2.1 Feedforward control without reflections 

Figure ^ displays a standard feedforward set-up for cancellation of random noise travelling at low- 
frequency in a duct (see Ref. Q). The duct is assumed to be rigid and either infinite or with anechoic 
terminations. A directional microphone detects the incoming noise (and not the secondary noise from 
the control source). An adaptive controller W filters this reference signal to generate appropriate input 
u at the control source to minimize noise at the error sensor. Below the first cut-off frequency of the 
duct, the noise e at the error sensor can be ideally written in the frequency domain as the sum of the 
primary and secondary noises: 

e = e-j'^'-^'x + ^e-j""/' =Fx + Hq (1) 



where c is the speed of sound, p the air density, S the cross-sectional area of the duct and q the 
loudspeaker volume velocity. If the loudspeaker dynamics is neglected so that q is directly the output 
of the control filter, the noise e at the error sensor will be cancelled if: 

9 <> 

q = Wx = -FH ^x=- —e-^^'^^-'^'x (2) 
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It appears that the volume velocity q which cancels noise is the causal filtering of the signal x if the 
reference sensor is located upstream of the loudspeaker (i.e. L > I). In practice however the secondary 
path H, which has to be inverted for perfect control, includes the dynamics of the loudspeaker and 
the delays inherent to anti-aliasing and reconstruction filters in a digital adaptive controller. The 
detection sensor must therefore be located further upstream in order to provide a reference signal 
with an advance larger than the overall delay in the secondary path. The purpose of this paper is to 
show that this classical result, which can be found in the textbooks on ANC (cf. Ref. fT\), does not 
apply to single-channel cases involving noise reflections nor to multi-channel cases even in free field: 
sometimes the detection sensor must provide a reference signal with an infinite advance for perfect 
control. 

2.2 Feedforward control with one reflection 

In this section it is supposed that the transfer function H between the control source and the error 
signal is the sum of direct propagation, with delay Xi and arbitrary magnitude 1, and of one echo 
with delay X2 (T2 > '^i) and complex magnitude a resulting from one (possibly non-rigid) acoustic 
reflexion or from the combination of several reflexions arriving simultaneously at the error sensor : 

H{(o) = e-j"'^'+ae-^'^^' (3) 

If the direct sound is louder than the echo (|a| < 1), the inverse of the transfer function can be 



expanded as: 
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Therefore, in the time domain, has an Infinite Impulse Response including a non-causal part 
with finite length Xi. As in the infinite duct case, the optimal control filter —FH^ will be fully causal 
if the transfer function F from the reference sensor to the error sensor includes a delay larger than 
the propagation time from the control source to the error sensor. When the primary noise is a pulse, 
the infinite time length of can be interpreted in a similar way to that applicable to cross-talk 
cancellation (cf. Ref. (T5\ ): the expansion in @ shows that the primary noise is cancelled by the 
direct sound resulting from a pulse in the control source, and that the echo of the control pulse is in 
turn cancelled by the direct sound of a second control pulse; and so on (the echo of the n^^ control 
pulse is cancelled by the direct sound of a («+ 1)"^ control pulse). Since the direct sound is louder 
than its echo, an infinite but converging series of pulses is required to achieve perfect cancellation. 
However the echo can be louder than the direct sound (|a| > 1). This may seem unlikely at first sight 
but a simple acoustic set-up with this feature is shown in section 12.31 The inverse transfer function 
can then be expanded as: 
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In this case it appears that has an infinite non-causal response, such that causal noise cancellation 
requires detection of the primary noise with an infinite time advance. With a pulse as primary noise, 
the theoretical non-causality of control can be interpreted from the expansion in the primary noise 
is cancelled by the echo of a counter-pulse in the control source, and the direct sound of this counter- 
pulse is pre-cancelled by the echo of a previous pulse, and so on. In the end an infinite converging 
series of control pulses in the past is required for noise cancellation of a single primary pulse. 
The simple analysis above shows that ANC is infinitely non-causal as soon as the echo is louder than 



the direct sound in the response. However equation ^ also demonstrates that, if the echo is very 
much louder, coefficient a is large and the series in ^ converges quickly. In this case the non-causal 
terms for large negative times may be neglected. Eventually control may be critically non-causal 
especially when the echo is slightly louder than the direct sound. 

2.3 Optimal non-causal control in a corner 

Figure |2| shows an ANC set-up where the secondary path includes an echo louder than the direct 
sound: two infinite rigid planes define in 3D a comer with a directional detection sensor, an omni- 
directional error sensor and a control source in the bisecting plane. With the notations from figure |2l 
and if the control source can be seen as a 3D monopole, which is the case of a loudspeaker at low 
frequency, the secondary transfer function H from the loudspeaker volume acceleration q to the error 
signal e can be computed by taking image sources into account: 
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with Li = L, L\ = (Li + Z)^ + and L3 = Li + 21. If L2 is less than 2L\ then the first echo in 
the secondary path is louder than the direct path. It can be shown that this condition amounts to 
L > l+fil. Fi gure |3] displays the secondary path as a function of frequency when c = 340m/s, 
L = Im and / = 0.8m (in fact these distances have been slightly adjusted in the subsequent numerical 
simulations so that noise propagation over L and / takes an integer number of samples). Figure 
|4] displays the corresponding exact Impulse Response and its digital approximation computed by 
Inverse Fast Fourier Transform for a sampling frequency of 16384//z and a frequency resolution of 

When the incoming primary noise is a plane wave propagating along the bisecting plane, as sketched 
in figure 121 the transfer function F from the directional reference sensor (detecting the incident plane 
wave only) to the omnidirectional error sensor (detecting the incident and reflected waves) is given 



for the primary noise by: 
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In this equation A denotes the distance between the reference sensor and the error sensor, as shown in 
figure 111 Without taking the dynamics of a real loudspeaker into account, perfect noise cancellation 
at the error sensor will be achieved if the input to the control source is equal to the reference signal 
filtered by the optimal filter W: 



With two echoes in the secondary path, no series expansion similar to that of equation ^ can easily 
express the causality of the optimal control filter W . However the Impulse Response of W can be 
computed numerically by IFFT. Figure |5l displays the Impulse Response of W when A = Om for a 
sampling frequency of 16384//z and a frequency resolution of -^Hz- In the case without reflections, 
A > would have enforced the control causality but, in the corner, the response includes a very long 
non-causal component, with significant coefficients at negative times much larger than any propaga- 
tion time between the elements of the set-up in figure |2l This means that noise cancellation could be 
achieved only by moving the detection sensor far away from the corner (i.e. by enforcing A >> 0), 
which may not be possible in practice (e.g. if the primary noise source is not far away). 

2.4 Optimal causally-constrained control in a corner 

In the previous section it has been shown that perfect noise cancellation can require very long non- 
causal filters. In practice however, the aim of ANC is only to provide significant noise attenuation. 
To this intent the performance of causal control must be determined as a function of the advance with 
which the reference signal is provided to the feedforward controller. 

As a first try, the quantitative effect of the causality constraint can be assessed by computing the noise 
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attenuation achieved with the causal control filter obtained discarding the non-causal coefficients in 
the time response of the optimal filter. For the set-up shown in figure|2lwith L = Im and I = 0.8m, 
figure |6l displays the noise attenuation achieved with this causally-truncated filter, when A = Om, 
A = 136m and A = 1700m, which amounts to providing the reference signal with an advance of 
respectively Os, 0.4s and 5s. Figure|6lshows that a huge advance in the reference signal is required to 
achieve substantial active noise reduction with a causally-truncated control filter. 
However the causal truncation of the non-causal optimal control filter is not the optimal causally- 
constrained control filter. The causality constraint does not restrain the control performance as much 
as figure |6l suggests. The optimal causal filter depends on the predictability of the primary noise. 
Following Ref. 1 16 1, it can be written in the case of white noise as: 



where H = HminHaii is the factorization of the secondary path FRF H into minimum-phase and all- 
pass components and //„//, and {}+ means that the non-causal component of the quantity en- 
closed by the braces is discarded. As in the previous section, F denotes the transfer function from the 
detection sensor to the error sensor for the primary noise and H the secondary path transfer function. 
Figure displays the noise attenuation achieved with the causally-constrained optimal control fil- 
ter when the reference signal is provided with an advance of Os, = 0.03125s and 0.4s, which 
amounts respectively to A = Om, A = 10.625m and A = 136m. For this figure the factorization of 
H into minimum-phase and all-pass components was performed in the frequency domain with the 
complex-cepstrum method (cf. Ref. lO)- For a given advance in the reference, the optimal causal 
filter performs much better than the causally truncated non-causal filters: with a 0.4s advance, op- 
timal causal control leads to a large noise reduction as in the case of causally-truncated non-causal 
control with a 5s advance. However the reference signal must still be provided with a very large 




(9) 



advance for efficient control. Figure [8l displays the Impulse Response of the causal filter for A = Om 
and A = 136m. For A = Om it shows that the optimal causal filter is very different from the causal 
truncation of the optimal non-causal controller, but for A = 136m it confirms that, with a very large 
advance in the reference, the response is close to the time-shifted response of the optimal non-causal 
controller (displayed in figure|5l). 

In practice ANC implementations very often rely on adaptive transverse control filters whose co- 
efficients are optimized recursively in the time-domain with the Filtered-reference(X) Least-Mean 
Square algorithm (cf. Ref. (7||). Figure|9ldisplays the noise attenuation achieved with a 1024-tap con- 
trol filter continuously updated with an FXLMS algorithm for A = Om, A = 10.625m and A = 21 .25m, 
and primary white noise. The propagation time over A = 21.25m corresponds to the whole length of 
the control filter impulse response. The FXLMS converged very slowly for these results. Even after 
selecting the best convergence coefficient, several seconds were typically required before achieving 
noise attenuation, and leakage (cf. Ref. (T6i ) did not improve convergence. For A = 10.625m the at- 
tenuation is of the same order of magnitude as for the optimal causal filter computed in the frequency 
domain. The zeroes in the attenuation are at the same frequency as the zeroes of the secondary path 
H or of the reference signal spectrum \F\^. At these frequencies the FXLMS converges with diffi- 
culties (cf. Ref. It could be expected that detecting the reference signal at a distance larger 
than A = 10.625m improve control performance, as for the causally-constrained filters derived in the 
frequency-domain. However the poor result shown in figure |9l for A = 21.25m means that the opti- 
mum advance of the reference signal does depend on the length of the filter response in the case of 
FXLMS control; the best control performance is not achieved with the largest pre-delay when the 
control filter impulse response is too short. 

Figure^! shows the impulse response of a FIR control filter obtained when the adaptation process is 
frozen after minimization of the error signal. Figure displays the noise attenuation theoretically 



achieved with this specific control filter. The computation is performed in the frequency domain for 
purely white primary noise. The attenuation for A = 10.625m is clearly worse than the attenuation 
displayed in figure|9lfor a continuously adapting control filter. Therefore the permanent adaptation of 
the control filter helps in ensuring good performance of FXLMS control, even after the residual error 
has apparently converged; a linear control FIR filter with permanent coefficient adaptation is in fact 
a non-Unear control filter which sometimes performs much better than any stationary linear control 
filter (see Ref. |^). 

2.5 Precautions against non-causality 

The series expansion in equation ^ suggests that control is critically non-causal when the echo in 
the secondary path response is only slightly louder than the direct sound. For the set-up of figure |2l 
the first echo is as loud as the direct sound if L = ^^/^ l ~ 1.2153/. With I = 0.8m this limit amounts 
to L 0.9722m, which is close to the value L = Im chosen for the numerical simulations. 
Figure Eldisplays the impulse response of the non-causal optimal filter W = —H~^F for L = 0.95m. 
The control filter has a very long response with only a relatively short (but non-zero) non-causal 
component. Figure displays the impulse response for L = 1 .05m. In this case the non-causal 
response decays more rapidly than for the case L = Im. These figures show that a small change in 
length L can dramatically affect the causality of optimal control. Control is highly non-causal only 
for the values of L which lead to a first echo slightly louder than the direct sound. 
In practice ANC set-ups often involve secondary paths with complicated time response including 
many echoes. It is not easy to detect without measurements if a given set-up gives rise to non-causality 
because of reflections, even if, as the previous simulations suggest, this will be critical only for very 
specific arrangements. One simple step towards preventing infinite non-causality in the single-channel 
case with reflections is to move the control source and the error sensor as close as possible to each 



other. In this way the direct sound is maximized. Furthermore reducing delays in the secondary path 
improves the convergence of adaptive algorithms such as the FXLMS (cf. Ref. |16|). However in 
practice other constraints may impose a minimum distance in the secondary path. For example the 
error sensor must usually be in the far-field of the control source to achieve noise reduction over 
more than a very narrow area. In the end infinite non-causality must be remembered as a theoretical 
possibility in active control of broadband noise. 

3 Non-causality in the multi-channel case in free field 
3.1 Optimal Control in a 2x2 case 

Figure shows a set-up for ANC in free field with two secondary sources and two omnidirectional 
minimization microphones. The incoming primary noise is assumed to be a plane wave. The loud- 
speakers are idealized as monopoles with volume accelerations qi and q2- The microphone #1 is 
taken as the reference sensor for feedforward control. The feedback from the control sources to the 
reference signal x will not be considered in the derivation of optimal control. In practice this can be 
done if the controller has an Internal Model Control structure including cancellation of the feedback 
path (see Ref. fl9l ). It is also assumed that the output control filters Wi, W2 with input x are directly 
the volume accelerations qi, q2 of the control sources, the dynamics of the loudspeakers is not taken 
into account. 

With the notations from figure (141 and with k = (n/c, the total noises p\, p2 at the minimization 
microphones can be written as: 
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control filters Wi, W2 are then given by: 
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It can be assumed, without loss of generality because the indexes can be switched, that d^d^ < d2dj,. 
After introducing a = did4/d2d3 and D = d2 + d^, — d\ — d^, equation (fTTl can be re- written and 
expanded as: 
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The causality of perfect noise cancellation can be discussed from this expression. Firstly, g-i^'^scose 
is a time delay if cos 6 > 0, that is if the primary noise is impinging the reference microphone before 
all the minimization microphones. In the opposite case ^-i^^'^scose (jg^otes a finite time advance equal 
to the propagation time of the primary noise between the reference microphone and the minimiza- 
tion microphones which receives it first. Secondly, the matrix in equation (I12t includes finite time 
advances. If D > one has di + d4 — d2 < dj, and di + d4 — dj, < d2, which means that these time 
advances are always smaller than the largest propagation time in the secondary path. In this case 
where D > 0, the series expansion in equation (fT2t will also be fully causal. Therefore, if D > 0, the 
single-channel result without reflections of section ITTI generalizes : noise cancellation is causal if the 
incoming primary noise is detected before it impinges the first minimization microphone, with a time 
advance larger than the largest propagation time in the secondary path. However, if D < 0, the series 
expansion in equation (fT2l gives rise to an infinitely-long non-causal impulse response for optimal 
control. This can occur indeed, for example in the case dy = 2.8m, d4 = 8.8m, and d2 = d3 = 5m. By 
inspecting the case did4 > d2dT„ it finally appears that noise cancellation will be infinitely non-causal 



in the general case if: 

{d\dn — d2d3){d\+ di^ — d2 — d-i) <Q (13) 

Figure[^displays four arrangements with increasing distances d\ and J4. At first sight it is not easy to 
detect significant differences between the arrangements but, according to equation (fT3l . the arrange- 
ments of figure EJa) and fTSl b) lead to infinitely non-causal noise cancellation. Figure fT6l displays 
the Impulse Responses of the optimal control filters, computed by IFFT with a sampling frequency 
of 16384Hz, when the primary noise is a plane wave with incidence angle = 0°. The response for 
case [TSla) is non-causal but it decays rapidly for negative times. Case ITsT b) leads to a very long 
non-causal response, which illustrates the non-causality of the series in equation The response 
for case[l5lc) has a finite short-length non-causal component but the causal component is very long, 
which would be a problem for practical implementation, whereas the response for case ITSl' d) should 
pose no problem in practice. 

The criterion given in equation ^Jcan be interpreted by considering active control in the 2x2 case as 
the combination of control with two individual single-channel devices and of cross-talk cancellation. 
Firstly, parameter a = did4/d2d^ can be seen as the averaged magnitude of the cross-talk at the 
minimization sensors. It quantifies the average ratio between the "direct" sound received by each of 
the sensor by the corresponding single-channel source and the "cross-talk" sound coming from the 
other source; if the arrangement of sensors and actuators is fully symmetric {i.e. di = and d2 = dj), 
a exactly measures this ratio at the two sensors. Secondly, parameter D/c = {d2 + d^ — d\ — dj^) / c 
measures the averaged difference between the times of arrival of the "direct" and "cross-talk" sounds 
at the minimization microphones. If a < 1 and D > 0, in average the "direct" sound is louder than the 
cross-talk and it reaches the minimization microphone before. If a < 1 and D <0, the cross-talk is not 
louder but it has to be cancelled before the direct sound, which requires an infinitely long non-causal 
response. If a > 1, the indexes of the sources can be switched for interpretation but, alternatively. 



the cross-talk can be seen as an echo louder than the direct sound if D > 0, and the discussion of 
section IT2I on control of an echo louder than the direct sound applies. Although there is no acoustic 
reflection in the free-field set-up of figure El the cross-talk can eventually be seen as an echo whose 
control implies non-causality if it is louder than the direct noise, as in the single-channel case with 
reflections. 

3.2 Causally-constrained control in a 2x2 case 

Figure^Jdisplays the noise attenuation achieved with causally-truncated control for the arrangement 
of figure fTSl b). For these simulations in the frequency domain, the reference signal was fed to the 
controller with the additional time advances i = d\/c = 8.24ms, x = Is and x = 2s. Figure (TTlshows 
that very large advances are required for large attenuations. The finite advance x = /c is the mini- 
mum advance required for the secondary noise to reach all the minimization microphones before the 
primary noise; it compensates the finite time advances which appear in equation (fT2l . If optimal con- 
trol were not infinitely non-causal, this finite time advance would be large enough for optimal control 
to be perfectly causal, and an infinite noise attenuation should be achieved. The very poor attenuation 
shown in figure [n] for x = d\/c illustrates the dramatic effect of non-causality upon control perfor- 
mances. In figure [n] the noise attenuation is higher at minimization microphone # 2 because of the 
propagation time of the primary noise between the two microphones, which is extra time for causal 
control. 

In the multi-channel stationary case, the optimal causally-constrained control filter can be written 
as in equation © as a function of the minimum-phase and all-pass components of the matrix of 
secondary paths llT6l . However, spectral factorization, which is an intermediate step in computing the 
optimal controller, is considerably more difficult to perform numerically in the multi-channel case 
than in the single-channel case (cf. ll20l ). Instead, in the context of active control, it is more efficient 



to compute sub-optimal controllers using adaptive algorithms in the time-domain. Fi sure fTSldisplavs 
the noise attenuation achieved with a 1024-tap control filter continuously updated with an FXLMS 
algorithm, for several additional advances x in the reference signal. The FXLMS took several seconds 
to converge to the results shown in figure^] For x = 8.24ms, control performs poorly and the noise 
is not reduced at all at microphone #1. For x = 31.25ms some noise reduction is achieved at the two 
minimization microphones. However, in this case, the noise reduction is low at some frequencies; it 
can be checked that the matrix of secondary paths is poorly conditioned at these frequencies, which 
is known as slowing down the convergence of the FXLMS algorithm (see Ref. [16|). Finally, the 
noise reduction is poor again for x = 62.5ms. As in section l2!4l for the single-channel case, this shows 
that, when optimal control is largely non-causal, the performance of Active Noise Control using a 
relatively short-length FIR filter depends strongly on the timing with which the incoming noise is 
detected. 

3.3 Precautions against non-causality 

In section 1231 it was shown that, in the single-channel case, non-causality was especially critical for 
arrangements where a slight change in the geometry would change the control filter with a non-causal 
response into a control filter with a long causal response. For the multichannel case of figure the 
series in equation (IT2l shows that the filter Impulse Response is long when a pa 1, which means that 
the secondary path H is close to a singular matrix at all frequencies. Therefore, in more complicated 
cases, the condition number of the secondary path matrix at low-frequency may be an indicator of 
possible non-causality. 

Figure (191 displavs the condition number at low-frequency (i.e. for A: — > 0) of the secondary path 
matrix H, as a function of parameter di when d4 = di+ 6m and d2 = d^= 5m. The vertical hnes in 
the plot delimit the values of d\ for which control is infinitely non-causal according to equation dTSl . 



The four cases shown in figures and are marked out by stars in figure It appears that very 
long non-causal Impulse Responses were obtained for the control filters when the condition number 
of the secondary matrix was high. In this 2x2 case, good conditioning appears to be a protection 
from the critical non-causality arising from the inversion of the secondary path matrix. Furthermore, 
for adaptive feedforward control, good conditioning of the secondary path is required for efficient 
convergence of the FXLMS algorithm (see Ref. (TEi ). In section 1231 a connection was also made in 
the single-channel case with reflections between causality and FXLMS convergence. From these two 
simple cases it may be inferred in practice that poor convergence of the FXLMS algorithm may be 
the sign of optimal control non-causality. 

As in the single-channel case, non-causality can be prevented in the case of figure by bringing 
the loudspeakers and the error microphones in closer pairs. In this way the condition number of the 
secondary path matrix is reduced because the diagonal of the matrix is made dominant. Regularization 
techniques can also be used to dampen the effect of poor matrix conditioning, e.g. by including a 
leakage coefficient in the FXLMS adaptation formula (see Ref. fT6ll ). Finally, the author believes that 
improving the condition number of the matrix of secondary paths also helps in preventing control non- 
causality in multi-channel set-ups involving more than two actuators and sensors. However, although 
non-causal set-ups with more channels can be built by adding sensor-actuator pairs to the set-up of 
figure fTSl' b'). a connection between causality and matrix condition number could not be made in the 
general multi-channel case. 

4 Conclusion 

In this paper it has been shown that active cancellation of broadband noise can theoretically require 
control filters that include an infinitely non-causal impulse response. In the single-channel case, 
this can happen because of high-level sound reflections. In the multi-channel case, the analysis of 



a 2x2 set-up has shown that this can occur even in free field because of the algebraic inversion of 
the secondary path matrix required for perfect cancellation. The cross-talk cancellation in the multi- 
channel case can be seen as equivalent to the cancellation of primary noise echoes in the single- 
channel case. Numerical simulations in the time and frequency domains have shown that, in practice, 
situations are possible where a significant noise reduction cannot be achieved without detecting the 
incoming primary noise with an advance very much larger than the propagation time delays in the 
secondary path. 

The inverse of the secondary path, which can have an infinite non-causal response, is required not 
only for Active Noise Cancellation but also for accurate reproduction of sound fields in applications 
such as "virtual acoustic imaging" (see Ref. |15|). In these applications, the computation of the 
inverse path need not be performed in real time, as is the case for ANC, and a processing time can be 
afforded which reduces the consequences of non-causality. However, designers of virtual acoustics 
devices must be aware that inverse acoustic responses can include time advances very much larger 
than the propagation times or the reverberation time in the direct response. 

It is not easy to prevent a large non-causality of the inverse acoustic responses in practical compli- 
cated cases, especially in rooms with many resonant modes. The analysis of the single-channel case 
in a comer and of the 2x2 case in free field have shown that one possibility is to mount actuators and 
sensors in close pairs, so that sound reflections and cross-talk components have lower levels than the 
direct sounds. This approach also improves the rate of convergence of adaptive control algorithms 
such as the Filtered-Reference Least Mean Squares. It also fits in with the current trend of designing 
ANC devices which are decentralized, in order to reduce the complexity of systems with many chan- 
nels. However the sensors must usually be in the far field of the actuators for global active control. 
Therefore infinite non-causality cannot always be prevented, and must be kept in mind as an ever 
possible limitation for active control of random noise. 
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List of figure captions 

Figure n A typical feedforward Active Noise Control set-up without reflections 
Figure 121 A feedforward Active Noise Control set-up in a comer 

Figure 0] Impulse Response of the secondary path, * continuous time response, — discrete-time 
model. 

Figure |3l Frequency Response Function of the secondary path. 
Figure 121 Impulse Response of the optimal non-causal control filter 

Figure |6l Noise attenuation with causal truncation of optimal control filters, — without additionnal 
advance in the reference signal, - - with a 0.4s advance, ... with a 5s advance. 

Figure Noise attenuation with causally-constrained optimal control filters computed in the fre- 
quency domain, — without additional advance in the reference signal, - - with a 31.25ms advance, ... 
with a 400ms advance. 

FigurejSl Impulse Response of causally-constrained optimal control filters computed in the frequency 
domain, ... without additional advance in the reference signal, — with a 0.4s advance. 
Figure 121 Noise attenuation with a 1024-coefficient control FIR filter updated with an FXLMS algo- 
rithm, — without additional advance in the reference signal, - - with a 31.25ms advance, ... with a 
62.5ms advance. 

Figure QOl Impulse response of a 1024-coefficient control FIR filter frozen after minimization of the 
error signal, with a 3 1 .25ms advance in the reference signal. 

Figure Noise attenuation with 1024-coefficient control FIR filters frozen after minimization of 
the error signal, — without additional advance in the reference signal, - - with a 31.25ms advance, ... 
with a 62.5ms advance. 

FigureEl Impulse Response of the optimal non-causal control filter for L=0.95m. 
Figure^l Impulse Response of the optimal non-causal control filter for L=1.05m. 



Figure [m A 2x2 ANC set-up in free field. 

Figure 4 ANC arrangements leading to optimal filters very different with respect to causality. 
FigureEl Impulse Responses of the optimal control filters for the 4 arrangements shown in figurefTSl 
Figure El Noise attenuation with causal truncation of optimal control filters, — with an 8.24ms 
additional time advance in the reference signal, - - with a Is advance, ... with a 2s advance. 
Figure Noise attenuation with a 1024-coefficient control FIR filter updated with an FXLMS 
algorithm, — with an 8.24ms additional time advance in the reference signal, - - with a 31.25ms 
advance, ... with a 62.5ms advance. 

Figure El Condition number of the secondary path matrix at low-frequency as a function of geomet- 
ric parameter di,ford2 = ds = 5m and d4 = di+ 6m. The stars mark out the values corresponding to 
the arrangements diplayed in figure (TSl 
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Secondary path in the frequency domain 
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